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ABSTRACT 

The  objective  of  this  work  has  been  to  study  the  interaction  between  heat  release  and 
mixing  in  compressible  shear  layers  by  analysis  and  computation,  with  an  eye  to  finding 
flow  configurations  that  maximize  the  heat  rdease  per  unit  distance  in  the  stream  direction. 
The  principal  prediction  is  that  heat  release  asymmetry  across  the  layer  can  enhance  mixing 
over  the  non-heat  release  case,  but  the  effect  appears  too  small  to  yidd  practical  benefits  at 
this  time.  Time-dependent,  three-dimensional  numerical  simulations  of  a  shear  layer  with 
weak,  steady  heat  release  have  shown  that  such  heat  release  need  not  decrease  mixing,  but 
the  expected  increases  are  also  absent,  or  too  small  to  be  detected  currently.  However,  the 
original  non-heat-release  theory  has  been  successfully  extended  to  predict  the  mixing 
behavior  of  three-dimensional  planar  layers  and  round  compressible  jets.  This  bolsters 
confidence  in  the  generality  of  die  principles  underlying  the  analysis.  Invoking  a  mixing 
maximum  principle,  the  extended  theory  gives  a  satisfectory  analytic  expression  for  mixing 
ratio  when  s  3.  A  variational  formulation  of  die  heat  rdease  problem  with  a  functional 
dependent  on  the  square  of  streamline  curvature  has  proved  intractable. 


NOMENCXATURE 

a^  s  speed  of  sound  in  reference  fiame  of  large  structures  of  layer 

a.  =  speed  of  sound  in  the  free  stream 

b  =  parameter  in  expression  for  heat  release,  hf 

c  =  parameter  in  expression  for  heat  release,  hf 

hf  =  enthalpy  due  to  heat  release  in  the  shear  layer 

m  =  mixing  rate 

m  =  non-dimensional  mixing  rate 

m^  =  reference  mixing  rate 

m^  =  reference  iiKompressible  mixing  rate 

=  large  structure  convection  velocity 
q_  =  reference  velocity 

s  =  refereiK:e  surface  area  for  mixing  rate  determination 

t  =  time 

u  =  streamwise  perturbation  velocity  due  to  shear  layer  s  +  fj, 

V  =  normal  perturbation  velocity  due  to  shear  layer  =  -  Fx 

x  s  streamwise  coordinate 


y  =  coordinate  transverse  to  the  mixing  layer 
D  =  parameter  in  expression  for  heat  release,  hf 

Cy  =  specific  heat  at  constant  volume 

Ey  =  perturbation  energy 

E,  =  energy  entrained  due  to  mbdr^ 

Eg  =  reference  energy 

E.  =  total  energy 

F  s  square  of  the  streamline  curvature, 

H  =  total  (volume  integrated)  heat  release  per  flow  distance 
K  s  streamline  curvature 

Me  =  convective  Mach  number  of  large  structures  in  a  shear  layer 

M.  «  reference  Mach  number 

NL  =  jet  exit  Mach  number 

R]^  s  rotational  kinetic  energy  of  a  vortex 

^  =  streamline  ardength 

T  =  average  layer  temperature 

T.  =  reference  temperature 

Uq  =  reference  velocity 

V  =  entrairunent  vdodty 

V  =  control  volume  in  mixing  analysis  wifii  no  heat  rdease 
a  =  parameter  in  odd  component  of  heat  release,  hf 

=  l-K2,orl-H? 
y  =  specific  heat  ratio 

1)  s  transverse  source  coordinate  in  Poisson  equation  integral  solution 

n'  =  n-y 

0g  -  hdicalar^e  of  a  helical  jet  vortex 

p  s  fluid  density 

$  =  potential  for  inotational  part  of  perturbation  vdodty 

yr  =  stream  fimetion  for  solenoidal  part  of  perturbation  vdodty 

u  =  vortidty 

£  =  streamwise  source  coordinate  in  Poisson  equation  integral  solution 

r  =  vortex  strength 

j 

INTRODUCTION 

A  major  problem  in  developing  hypersonic  aircraft  is  the  design  of  an  efficient 
propulsion  system.  One  of  the  most  difiicult  aspects  of  propulsion  design  is  to  get  the  fuel 
and  air  mixed  suffidenfly  so  that  combustion  is  maximized  inside  the  engine.  At  the  mass 
flux  rates  of  fuel  and  air  necessary  to  produce  the  design  thrust,  the  fluid  streams  are 
moving  at  high  supersonic  speeds.  As  the  differer^e  between  the  stream  Mach  numbers 
increases,  the  mixing  of  the  streams  decreases,  with  a  corresponding  loss  in  combustion 
eflidency.  This  mixing  decrease  correlates  wdl  with  the  "convective  ^^ch  number,"  which 
is  essentially  a  density  weighted  difference  of  the  Mach  numbers  of  the  streams.  If  the 
convective  Mach  number  is  imity,  then  experimental  results^  show  that  the  mixing  is  about 
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25%  of  its  value  at  a  convective  Mach  number  of  zero.  There  is  a  considerable  body  of 
literature  on  the  causes  of  this  deterioration  in  mixing,  including  experimental, 
computational,  and  theoretical  studies. 

Earlier  work^  examined  the  mixing  of  supersonic  ^ear  layers,  both  theoretically  and 
computationally,  with  the  condtision  that  mixing  is  inhibited  at  higher  convective  Mach 
numbers  (M^  because  energy  that  would  be  directly  transmitted  into  mixing  at  M^  s  0 
instead  goes  into  nonmixing  gas  compression.  This  compressibility  effect  is  l^gely  two- 
dimensional;  introducing  three-dimensionality  into  the  flow  (streamwise  vortices)  r^uces 
but  does  not  eliminate  it  for  M^  >  .7.  Comparison  of  experimental  growth  rates  to  the 
theory  in  Reference  2  suggests  that  such  changes  in  flow  topology  are  occurring  naturally 
for  high  M^  layers,  and  little  else  can  be  done.  Without  these  changes  mixing  would  be  even 
less  than  is  measured. 

Dismaying  as  these  results  are,  it  should  be  borne  in  mind  that  the  ultimate  goal  is  not 
just  mixing  of  inert  species,  but  rather  mixing  of  fuel  and  oxidant  for  combustion.  The  heat 
release  associated  with  combustion  can  be  expected  to  modify  the  structure  of  the  shear 
layer,  and  this  in  turn  affects  its  mixing  rate.  Thus  it  is  appropriate  to  study  the  coupling  of 
mixing  and  heat  rdease,  for  there  may  be  interactions  between  them  that  could  lead  to  more 
effkdent  combustion  than  inert  mixing  results  alone  indicate. 

There  have  been  a  number  of  experimental  studies,^'^  as  well  as  some  niunerical 
work,^'^  on  the  effect  of  heat  release  on  low  speed  mixing  layers  (Mj^  <  .35).  All  studies 
conclude  that  heat  release  decreases  the  layer  growth  rate,  except  for  &e  work  of  Kdler  and 
Daily.  This  one  experiment  is  distinguished  by  its  use  of  hot  combustion  products  in  one 
stream  to  ignite  completely  premixed  fuel/air  in  the  other  stream.  Other  studies  assume  a 
symmetric  placement  of  unbumt  fuel  in  one  streanv  oxidant  in  the  oflier.  Kdler  and  Daly 
are  the  only  ones  to  report  an  increase  in  mixing  layer  growth  rate  with  increased  heat 
release.  Within  the  framework  of  the  present  analyses  such  asymmetries  are  demonstrated 
to  be  mixing  enhancers. 

The  previously  developed  theory^  for  nonreacting  shear  layers  produced  an  analytical 
prediction  of  the  reduced  mixing  observed  in  experiment.  This  prediction  was  made  by 
discarding  previous  arguments  relating  reduc^  mixing  to  reduced  instability^*^  and 
employing  the  insight  that  this  phenomenon  is  explainable  by  considering  wholly 
irrotational  compressibility  effects.  In  the  present  work  this  analysis  is  extended  to  include 
the  effects  of  weak  heat  release  and  entropy  generation  with  and  without  vorticity.  The 
theory  is  meant  to  be  valid  at  hi^h  speeds  -  0(1)),  and  this  results  in  somewhat  different 
governing  equations  than  appear  in  analysis  of  heat  release  in  flows  as  M^  0.’  The  results 
indicate  that  mixing  rates  can  be  either  increased  or  decreased  by  proper  tailoring  of  the 
heat  release  profile  across  the  layer. 


ELEMENTS  OF  THE  THEORY  FOR  INERT  MIXING 

The  physical  configuration  of  a  compressible  free  shear  layer  is  shown  in  Figure  1, 
where  6)2  vorticity  contours  have  been  drawn  from  the  results  of  twoidimensional  direct 
numerical  simulations  of  this  flow.  Two  parallel  streams,  initially  separated  by  a  splitter 
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plate;,  and  semi-infinite  in  directions  transverse  to  the  mean  flow,  mix  at  their  interfile.  The 
upper  stream  is  characterized  by  Pj,  U],  a^,  and  the  lower  by  P2,  U2,  a2  (U^  a  U2  assumed). 
Static  pressures  in  the  stream  are  matcl^  at  the  tndling  edge  of  the  splitter  plate  so  that  no 
expansion/con^ression  waves  contaminate  the  streamwise  devdopment  of  the  flow.  Initial 
instability  of  this  shear  layer  causes  it  to  roll  up  into  spanwise  vortices  which  subsequently 
pair  due  to  a  further  spatial  subharmonic  instability.  As  the  roll-up  and  growth  of  vortices 
proceeds,  the  width  of  the  shear  layer  grov^  downstream,  and  increasing  amounts  of  fluid 
are  entrained  from  both  sides.  The  growth  rate  of  incompressible  shear  layers  depmds  on 
their  velocity  and  density  ratios.  High  speed  shear  layers,  however,  grow  much  more 
slowly  than  their  incompressible  counterparts  at  the  same  velocity  and  density  ratios. 
Bogdanoff,^  and  later  Papamoschou  and  Roshko,^^  suggested  that  the  effect  of 
compressibility  on  growth  of  a  shear  layer  is  characterized  by  Mach  numbers  ^ 
measured  rdative  to  the  convective  speed  of  the  large  structures. 

Thus 


A  relation  between  and  is  obtained  by  assuming  that  the  two  streams  have 
equal  static  pressures,  and  that  their  stagnation  values  in  the  reference  frame  are  also 
equal.  If  it  is  further  assumed  that  the  ratio  of  spedfic  heats  y|  for  each  stream  are  equal, 
then  Me  s  M^^  -  M^,  and  the  convective  velocity  may  be  c^culated  as  a  sound  speed 
weightra  average 


Vi  *  Va 
*1  ♦  *2 


(2) 


The  dominant  feature  of  a  mixing  layer  at  lower  values  of  M^  is  the  existence  of  large 
vortical  structures.  If  these  structures  are  regarded  as  characteristics  of  the  flow,  then  an 
analysis  based  on  the  behavior  of  these  structures  wifli  increasing  convective  Mach  number 
should  also  reflect  the  behavior  of  the  complete  mixing  layer  with  increasing  M^.  The 
theory  of  Reference  2  has  three  elements.  The  first  is  an  expression  for  the  energy  wifldn  a 
control  volume  circumscribing  a  portion  of  the  turbulent  shear  layer,  in  terms  of  flie  mixing 
rate  (rate  of  increase  of  mass  within  the  volume)  and  convective  Mach  number  M^.  In  the 
second  element  the  mixing  rate  of  an  isolated  vortex  is  shown  to  be  proportional  to  the 
square  root  of  its  energy,  and  this  establishes  a  scaling  between  the  mixing  rate  due  to  an 
arbitrary  vortex  and  that  of  a  reference  vortex,  based  on  the  ratio  of  their  energies. 
Identifying  the  energy  and  mixing  rate  of  the  control  volume  analysis  with  that  of  an 
arbitral  vortex  yields  a  relation  between  the  nondimensionalized  mixing  in  the  shear  layer 
and  the  energy  and  mixing  due  to  a  reference,  isolated  vortex  sitting  in  a  uniform  tree 
stream  of  Mach  number  M^  The  basic  result  can  be  written  as 
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where  m  denotes  non-dimensional  mixing  of  the  mixing  layer,  the  subscript  i  indicates 
values  at  and  mo  doiotes  non-dimensional  mixing  associated  with  a  refoence  vortex 

in  which  the  rotational  Idnetk  energy  GUCE)  is  GiKElo- 

Finally,  in  the  third  element  of  the  theory  the  reference  vortex  is  taken  to  be  the  bound 
circulation  of  a  lifting  airfoil  Both  the  mixing  and  energy  of  this  reference  "vortex"  can  be 
calculated  analytically  from  solution  of  the  Prandd-Glaueit  equation  with  M.  =  M^. 

=  0  (4) 

The  analytic  pre^ction  of  the  theory  for  the,mixing  ratio  at  lower  values  of  M<.  (M^  s  0.7) 
(ratio  of  mixing  m  at  to  mixing  at  Mc=0,  mi,  holding  density  and  vdocity  ratio  constant) 
is 


(1  - 


(5) 


Plotted  as  the  solid  line,  the  central  character  of  this  relation  with  respect  to  the 
experimental  scatter  for  <  .7  can  be  seen  in  Figure  2.  A  modification  to  this  theory, 
described  in  later  sections,  extends  it  to  greater  values  of  M^,  by  allowing  for  three- 
dimensionality 

The  airfoil  analogy  and  its  solution  by  potential  methods  is  a  framework  that  embodies 
several  important  physical  efrects  relevant  to  mixing.  Replacement  of  the  shear  layer  by  the 
velocity  difference  across  the  airfoil  reduces  the  mixing  problem  to  the  simplest  terms: 
given  a  horizontal  velocity  difference  (whose  integral  over  the  chord  is  the  bound 
circulation  of  the  airfoil),  what  is  the  vertical  velocity  (downwash,  entrainment  velocity)  that 
it  induces?  The  answer  to  this  question  provided  by  a  potential  flow  solution  to  the  airfoil 
problem  necessarily  implies  that  the  predicted  effects  are  irrotational  The  quantitative 
accuracy  of  Equation  (2)  in  predicting  observed  mixing  decreases  strongly  implies  that  both 
the  instability^'*  and  vorticity  redistribution**  arguments  simply  miss  the  point.  The 
physical  explanation  for  the  irrotational  decrease  of  mixing  is  simpler:  compressibility 
directly  attenuates  the  large  structures'  ability  to  deflect  fluid  away  from  the  free  streams 
and  into  the  mixing  layer.  In  incompressible  layers  "pushing"  on  the  free  stream  by  the 
entrainment  velocity  deflects  fluid  into  the  lay^;  in  compressible  layers  "pushing"  produces 
fluid  compression  also  and  thus  less  deflection. 
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MIXING  THEORY  WITH  HEAT  RELEASE  AND  VORTICITY 


The  same  three  elements  make  up  the  theory  for  mixing  with  heat  release  as  did  the 
theory  for  mixing  without  heat  However,  the  simple  scaling  rdations  that  lead  to  Equation 
(5)  can  no  longer  be  obtained  v^en  heat  release  is  included.  Instead,  the  analysis  results  in  a 
general  expression  for  the  nondimensional  mixing  that  can  be  made  specific  by  assumption 
of  a  particular  he?'  release  configuration. 

The  analysis  is  best  understood  in  the  cont&ct  of  the  control  volume  and  vortex  model 
depicted  in  Figures  3  and  4.  The  mixing  per  unit  span  of  the  vortex  is  defiled  by 

m  =  I  p^lvldx  (6) 

where  now  the  entrainment  or  downwash  velocity  v  is  assumed  to  have  three  components 

v  =  Vpg  +  Vh  +  V4,  (7) 

This  decomposition  looks  ahead  to  the  inhomogeneous  Prandtl-Glauert  equation  that  will 
ultimately  be  solved  for  v,  and  which  has  homogeneous  solution  Vpg  and  specific  solutions 
vi,andv^  Then  the  mixing  can  be  written 

m  =  mpg+n^  +  m^  (8) 

where  is  the  mixing  due  to  heat  release,  m^  is  the  mixing  due  to  vortid^,  and  ntpg  is  the 

value  using  the  Prandtl-Glauert  equation. 

Similar  to  the  original  analysis,^  the  control  volume  balance  yields 


A  ^ 

where  Eh  is  the  eneigy  in  the  mixed  region  due  to  heat  release  and  Eg  is  the  total  energy  in 
the  mix^  region.  One  may  also  write  for  the  isolated  vortex 


m 


-  V 


1/2 


(10) 


A 

where  E^  and  (Eg  -  Eh)  are  the  rotational  kinetic  energy  of  the  Prandtl-Glauert  vortex  at  the 
reference  condition  and  the  current  condition  respectively.  Using  Equation  (9)  and 
proceeding  as  in  Reference  2  gives  the  equation 
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A  A  _  A  *1/0 

-  ®h'  “o  -  *h'  -  - 

772— 


E 


B 


where  E_aiKl  are  givett  by 


A 

B 


B 


Vo"" 


B  =  C  T  (  1  +  Illlil  i|2  ) 

a>  V  2  *0 


and  V  is  a  normalizing  volume  derived  from  the  control  volume  balance.  See  Reference  2 
for  details. 


Equation  (11)  la  quadratic  In 


thUB 


^®a  - 


1/2 
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1e  1/2 


± 


(12) 


and,  hence,  from  Equations  (9)  and  (11) 


m 


-2 

m 


4  r 

Hf 


i-f-**%* 


(13) 


o  o  o 

In  order  to  recover  the  nonreactLtg  result,  the  positive  sign  is  taken.  If  the  assumption  is 
made  that 


'“h 


+  m  I 
61 


<  < 


-  2 
m 
o 
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then  Equation  (14)  gives  the  simple  equation 


m  = 


o 

E 


4-  2  m  4  2m 


a 


6) 


(14) 


(15) 
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The  inert  mixing  theory  produces  just  the  first  term  on  the  right  hand  side  of  Equation  (15), 
and  this  leads  immediatdy  to  Equation  (3). 

To  proceed  further  the  entrainment  velocity  v  must  be  solved  for.  This  can  be  obtained 
from  a  Poisson  equation  for  the  irrotational  part  of  the  vdodty  field.  The  analysis  leading  to 
this  equation  is  quasi-steady.  While  molecular  mixing  within  the  layer  cannot  be  treated  in 
this  manner,  the  entrainment  process  by  which  the  layer  grows  can.  In  particular  note  that 
the  temporal  rate  of  change  of  large  structures  within  the  layer  is  of  the  order  of  the 
spreading  rate,  and  this  is  small.  The  steady  energy  equation  for  specific  enthalpy,  h,  of 
fluid  in  the  layer  is  given  by 


h+(l/2)<]P  =  lV,  +  h,  (16) 

where  h^  is  the  stagnation  enthalpy,  h|  is  the  heat  release,  and  q  is  the  magnitude  of  the  total 
fluid  velocity. 

If  the  inviscid  momentiun  and  energy  equations  are  combined,  the  steady  Crocco 
equation  results 


V(h  +  -^)  =  Vh<,  + Vh|  =  TVS  +  cpcw  (17) 

where  q  is  the  velocity  vector,  T  is  temperature,  S  is  entropy,  and  o  is  vortidty.  But,  h^  is 
taken  to  be  a  constant  and  thus 


Vhf  =  TVS  +  qxw  (18) 

In  the  absence  of  any  other  generators  of  entropy  the  heat  release  appears  as  an  increase  of 
entropy  and  vortidity.  Now  let  the  cartesian  p^turbation  velocity  components  u  and  v  be 
written  as 


u  =  ^x  +  Py;  v  =  (19) 

where  ^  is  the  perturbation  velocity  potential  and  yr  is  the  stream  funcbon  given  by 

rxx+ryy=-«/u,  (20) 

From  the  energy  equation  the  density,  p,  can  be  written  in  terms  of  the  velocity  and 
entropy  as 


'a  U  ^ 

OP  CD 

The  assumption  of  weak  heat  release  and  small  entropy  change 


/R)  (21) 
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(y-i)hj 


s  -  s 


«  1 


aUows  this  state  equation  to  be  linearized.  When  combined  with  the  continuity  equation 

V*(#>q)sO  (22) 

and  eliminating  second  order  terms  in  the  perturbations  (velocity,  heat  and  entropy),  a 
Poisson  equation  for  the  perturbation  potential  results. 


(1  -  M^)#  +  # 

-  XX  yy 


+ 

m 

•  Y 


The  right  hand  side  of  Equation  (23)  can  be  simplified  by  considering  the  two  components 
of  the  Crocco  equation.  At  this  order 

(hf),  =  TJ5, 

(24) 

(hf)y  =  T^-VjO 

Integrating  the  first  of  this  pair,  differentiating  by  y  and  comparing  with  the  second  shows 
that 


Then  defining 


a  T  S  -  I  U  « 

«P  J  ^ 


r  =  yi^  j  («/u^)  dy 


the  governing  equations  become 


(1  -  ir)$  +  #  = 

»  XX  yy 


(K-l) (h,). 


+  +  yM^r 

•  xy  '  «*  X 


and  Equation  (20).  The  formal  solution  of  the  above  equation  is  described  in  the  following 
section. 
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SOLUTION  FOR  THE  ENTRAINMENT  VELOCITY 


The  formal  solution  of  the  Poisson  Equation  (27)  for  the  perturbation  potential  can  be 
fouiui  in  any  text  on  mathematical  physics.  It  involves  integration  of  the  right  lumd  side  of 
Equation  (2^,  multiplied  by  the  appropriate  Great's  function,  over  the  spatial  domain  of  the 
problem.  In  the  context  of  the  norueacting  analysis  this  two-dimensional  volume  integral 
was  transformed  into  a  line  integral  over  the  chord  of  a  thin  lifting  airfoil  whose  bound 
circulation  was  given  by  the  velocity  difference  of  the  two  streams  of  the  mixing  layer.  Here 
the  same  approach  is  t^en,  with  die  "source"  function  on  the  right  hand  side  of  Equation 
(27)  assumed  to  be  confined  to  the  plane  y=0.  Using  the  integral  equation  analysis  of 
Reference  12,  a  formal  solution  to  Equation  (27)  can  be  written  as 

V  =  v„  +  Vi  +  V.  (28) 

where  v-^  denotes  the  value  for  the  Prandtl  Glauert  equation  (Equation  (27)  with  the  right 
hand  sicte  equated  to  zero).  V|,  and  v^  are  given  by 


o 


where 

L()  is  the  operator  on  If  given  by 


o 


and 


(30) 


(31) 


The  kernel  K  is  the  two  dimensional  Green's  function 
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K  =  -^lnl(x-{)2+i,2)l 


(32) 


for  the  upper  half  plane  D^,  and 


^=M£(ryy-y«/uj 


(33) 


The  structure  of  the  integrals  in  Equation  (31)  makes  clear  the  possibilities  for  mixing 
increase  or  decrease  when  heat  release  is  present  Specifically,  if  h^  is  even  or  >l  is  odd  in  y, 
then  Vj,  and  are  zero  and  the  mixing  is  unafiectecL  On  the  other  hand,  if  hf  is  odd  and  2  is 
even  then  both  additional  velocities  will  be  non-zero,  and  mixing  can  ^  increased  or 
decreased.  Vj,  can  be  either  positive  or  negative  dependent  upon  which  side  of  the  layer  the 
heat  release  is  concentrated,  while  the  sign  of  v^  depends  upon  a  combination  of  streamline 
curvature  and  vorttdty  whose  physical  interpretation  is  difficult  to  determine. 

Some  further  simplifications  of  the  expressions  for  entrainment  velocity  due  to  heat 
release  can  be  made  if  a  particular  mixing  modd,  the  vortex,  is  assumed.  In  Equation  (29), 
an  airfoil  analogy  is  used  and  the  chord  of  the  airfoil  is  assumed  to  be  unity.  If  the 
dominant  flow  structure  is  a  simple  vortex  centered  at  x^  then  this  chord  approadres  zero, 
or,  in  the  case  of  the  line  integral  in  Equation  (29) 

lxl>>UI  (34) 


Hence 


1.  (X)  =  I-  (x^) 

t.,  t,  o 


(35) 


Also,  then 

I  ^ 

0 

and  the  vdodty  V|,  is  given  by 


0  0 


(36) 


-"fi'V  ' 


(37) 


The  mixing  rate  m)|  is  then  given  by 
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s.  ■  ..-hH  I M  -M-  •  I  IH^r-M«l 


I  \ 


(38) 


(x)dx 


Ifl  can  be  rewritten,  using  an  integration  by  parts,  as 

*  5i  n  '  h^{{,-8)ldfd8 


(39) 


NUMERICAL  SIMULATIONS  WITH  STEADY  HEAT  RELEASE 


After  addition  of  an  independently  specifiable  heat  rdease  term  to  the  energy  equation 
in  our  Navier-Stokes  simulations,  the  eff^  of  steady,  but  spatially  varying,  heating  on  a 
two-dimensional  shear  layer  has  been  investigated.  The  heating  profile  is  specified  as  a 
shifted  Gaussian: 


hj(x,y) 


(40) 


By  var3ring  D,  b,  or  and  c  the  intensity,  displacement  and  spreading  rate  of  the  heating 
profile  can  be  adjusted  in  a  parameter  study.  Becaiise  of  the  particular  form  chosen  the  total 
heat  release  per  unit  distance  in  the  flow  direction  is  constant,  consistent  with  linear  growth 
of  a  reacting  shear  layer.  Previously  reported  experiment^  and  simulation^  have  left  the 
general  impression  that  heating  reduces  shear  layer  growth,  alfiiough  there  is  one  report^  to 
the  contrary.  Although  mixing  rate  increases  due  to  heating  asymmetry  were  not  found  in 
the  parameter  study,  neither  were  decreases  in  mixing  detected.  Figures  3  and  4  are  typical 
of  the  simulation  results  and  display  the  spreading  rates  and  temperature  contours  of  an 
Mc=.4  layer  with  and  without  heating.  Normalized  temperature  contours  ( (T-TJ/T. )  are 
shown  only  for  the  heated  layer,  since  at  the  same  contour  resolution  the  unheated  layer  has 
no  perceptible  variation.  The  parameter  values  for  the  heat  release  are  D=.002,  b=l.,  a=2. 
and  CS.OZ.  For  these  conditions  the  heat  is  released  completely  within  the  layer  and  the  time 
averaged  temperature  rise  of  the  fluid  within  the  turbulent  part  of  the  flow  reaches  a 
centerline  peak  more  than  30%  greater  than  its  unheated  value.  Yet  the  spreading  rate, 
measured  by  the  vortidty  tiiickness  is  essentially  unaffected.  Such  cases  have  not  been 
analyzed  in  detail,  but  they  indicate,  along  with  the  strong  heat  release  cases  reported  in 
Refd^nce  4,  that  there  can  heat  release  scenaric^  without  decreased  mbdi^. 
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AN  ANALYTIC  EXAMPLE  OF  HEAT  RELEASE 


Prompted  by  the  anomalous  results  of  the  simulations  an  analytic  heat  release  example 
was  solved  to  better  understand  the  functional  dependence  of  the  entrainment  velocity 
upon  heat  release  parameters.  The  field  equation  for  the  perturbation  velocity  potential  due 
to  heat  release  hfCx,y)  is 


2 

a 

OP 


(x.y) 

X 


(41) 


The  analysis  begins  with  the  expression  for  W,  tl^  vdodty  normal  to  the  layer,  gotten  from  a 
y-differentiation  of  the  expression  for  the  potential  $  that  solves  the  governing  Poisson 
equation.  This  is 


V(x,y)  =  ♦y(x,y) 


(y-1) 

2Jia^ 


f  1 

«  n 


_ ii. 


(x-{)^  + 


{y-u)^l 


di}  d(  (42) 


Note  the  asymmetry  of  the  kernel  in  Eq.  uidess  hf^  has  some  asymmetry  in  q  with 
respect  to  the  y-value  of  the  field  point,  the  q  integral  equates  to  zero.  Thus,  for  a  given 
heat  release  distribution  hf .,  it  is  only  its  component  odd  with  respect  to  y  tiiat  is  important 
to  the  entrainment  velocity  V  at  y. 


For  simplicity,  assume  a  separable  form  for  h{((,q) 


h^(^.q)  s  Hf({)g(q)  (43) 

where  (((),  g(q),  H  >  0  and  JJ  hf  dqd$  must  be  finite. 

Then 


V(x,y) 


(y-l)H 

2na^^ 


(y-q)g(q) 


[(x-{)^+  (y-q)  ] 


dq  d( 


q=- 


(44) 


Now  g(q)  can  be  decomposed  about  any  y  into  a  sum  of  odd  and  even  functions  go(y~Q)  + 
gg(y-q).  Only  g,,  contributes  to  the  integral,  and  so  Eq.  (44)  may  be  rewritten  replacing  g(q) 
by  sJy-tO-  To  understand  the  subsequent  behavior  of  V  we  choose  a  generic  form  for  gjy- 
il)  that  has  properties  appropriate  for  the  odd  component  of  a  single  maximum  (in  q)  heat 
rdeaseproffle.  It  is 
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Q^iy-n) 


-2a  (y-ip 

**  [1  +  (y-n)^  ^ 


(45) 


This  form  has  been  chosen  for  its  integration  properties,  its  structural  resemblance  to  the 
odd  component  (about  arbitrary  y)  of  the  fun^on  g(i|)  =  l/d+b^i}^),  and  for  its  property 
that  as  b^,  •*  -aS'(y~ii).  The  constant  a  controls  the  sign  of  the  contribution  of  the 

heat  release  above  and  bdo  w  y.  Substituting  Eq.  (45)  into  Eq.  (44)  yidds 


V(x,y). 


_ (y-»i)^ _ 

[(x-{)^+(y-i|)^]  [l  +  b^(y-ii)^] 


dq  (46) 


»1=— 

After  a  change  of  variable  to  z  =  y-i),  and  splitting  the  integrand  into  factors,  the  inner 
integral  can  be  evaluated  to  yield 


V(x,y) 


fl4>^(y-l)H 

^  2-2 
2xa  8 


[b(x-^)+l] 


(47) 


For  f^  0,  i.e.,  heat  release  increasing  in  the  streamwise  direction,  the  direction  of  the 
entrainment  velocity  is  controlled  by  the  sign  of  a.  This  in  turn  characterizes  whedier  the 
odd  component  of  the  heat  release,  at  a  given  y,  is  negative  above  or  bdow  y.  In  either  case, 
the  entrainment  velocity  goes  from  the  positive  portion  of  tiie  odd  component  towards  the 
negative.  This  situation  is  depicted  in  Figure  5. 

What  is  the  general  shape  of  go  for  a  single  maximum  heat  release  profile?  One  can 
construct  the  answer  graphically,  but  for 


9(q) 


1 

1  + 


2 

n 


it  can  be  demonstrated  that  for  q'  =  i^-y 


0^(0'  +  y)= 


_ -2n’Y _ 

[i  +  (ii*  +  y)^][i  +  <“*1*  +  y)^] 


(48) 


Thus,  for  positive  y,  the  negative  component  of  go  is  above  y,  and  the  oitrainment  velocity 
V  is  in  that  direction,  away  from  the  heat  release.  For  y  <  0  the  negative  portion  is  below  y, 
and  once  again  the  entrainment  velocity  is  away  from  the  heat  release,  llius,  for  increasing 
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heat  release  downstream,  >  0,  entrainment  is  decreased,  but  for  f^  <  0,  decreasing  heat 
rdeasedownstream>  entrainment  is  increased  by  heat  release. 


The  problem  with  the  numerical  simulations  now  becomes  clearer.  Despite  the  foct 
that  the  assumed  form  of  the  heat  release,  Eq.  (40),  is  not  separable,  the  streamwise 
derivative  (f^  can  be  calculated,  and  is  fcnind  to  vaiy  in  sign  with  bothx  and  y.  However,  it 
appears  to  be  uniformly  negative  for  regions  nearby  the  shear  layer,  which  means  the 
entrainment  velocities  are  towards  the  mixing.  The  magnitude  of  V  expected  in  the 
simulations  can  be  estimated  from  the  aiudysis  if  one  identiftes  Di^V(y-l)  in  (40)  with  H 
in  Eq.  (43),  with  b  in  Eq.  (45),  and  realizes  from  x-difrerentiatk>n  of  Eq.  (40)  that 


1 

(1  eae)^ 


is  the  maximum  value  of  the  streamwise  derivative  for  field  points  y  on  the  edges  of  the 
layer.  Thus,  hf^ "  [-1,  --03]  over  the  length  of  the  layer  and  we  estimate 


V 


-  o  (y-1) 

2  2 
2  na  p 

-  ^ 

2nl? 


Da 

m 


2 


(y-1) 


(2) (.002) 
“HTTSST 


-  -  .00076  (49) 

The  heat  release  distribution  induces  entrainment  vdodties  one  to  two  orders  of  magnitude 
smaller  than  the  typical  turbulent  values  of  V'  in  the  shear  layer.  Such  effects,  though  real, 
could  easily  be  in^tinguishable. 


EXTENSIONS  TO  INERT  MIXING  THEORY 

Despite  the  apparent  smallness  of  the  heat  rdease  induced  entrainment  vdodties,  the 
effects  are  believed  to  tbe  real  The  broad  validity  of  the  prindples  underlying  the  analysis 
can  be  demonstrated  by  extending  the  inert  mixing  theory  to  three-dimensions  and 
comparing  the  predictions  with  experimental  data.  The  way  in  which  three-dimensional 
effects  can  be  incorporated  in  the  theory  for  planer  layers  and  circular  jets  is  described  in  the 
following  paragraphs. 


PLANAR  SHEAR  LAYERS 

The  theory  developed  in  Ref.  2  is  based  on  two  main  argvmnents:  first  that  the  energy 
into  an  element  of  the  mixing  layer,  A6CD  in  Figure  6,  is  conserved;  secondly,  that  a 
percentage  (100%  in  two  dimensions)  of  this  energy  is  converted  into  rotational  energy 
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associated  with  the  dominant  mixing  mechanism,  which  is  assumed  to  be  a  two- 
dimensional  vortex.  This  vortex  is  represented  by  the  potential  boimd  vortex  associated 
with  a  lifting  wing  and  leads  to  a  simple  analytic  re^t  if  the  potential  vortex  can  be 
represented  ^  the  two-dimensional  PrandtlGlauert  eqpiation 

^#*x  +  #yy  =  0  (50) 

wl^re  #  is  a  perturbation  velocity  potential  aitd 

(51) 

The  basic  equations  derived  in  Ref.  2  for  the  nondimensionalized  mixing  rate,  m, 
defined  by 


m  =  m/p,,UoS  (52) 

where  p^,  Uq  are  the  reference  values  of  density  and  velocity  respectively  associated  with 
and  is  the  mixing  role  through  a  surface  area  s,  are  as  follows:  the  energy  entrained 

in  time  t,  is  given  by 


E.  =  C,Tm[l  + 


(53) 


where  Cy  is  the  specific  heat  at  constant  volume,  and  T  is  the  average  temperature^  the 
mixing  layer;  thto  is  Eq.  (9)  of  Ref.  ^  with  the  nondimensionalized  ntixing  rate,  m,  and 
instead  of  m  tiie  energy  of  the  vortex  E  is  given  by  Eq.  (13)  of  Ref.  2,  which  is 


(54) 


where  the  subscript  ”o"  denotes  a  known  reference  state.  Provided  M  is  such  that 

M^<<2/(y-l)  (55) 

it  is  shown  tit  Ref.  2  that  over  the  control  volume  V, 


A 

E 
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(56) 


and 
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(y-1) 
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q^dv 


j  v-p  * 

V 


(57) 
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where  M  in  Eq.  (55)  is  the  oncoming  Mach  number  of  a  t%vo-dimensional  vortex  and  in 
Eq.  (56)  is  the  square  of  the  perturbation  velocity  due  to  the  vortex,  superimposed  as  a 
heestream  of  energy,  E^  and  a  velocity  Eq.  (53),  (54),  (56),  and  (57)  are  genei^  relations 

and  are  not  restricted  to  tjyo  dimensions.  In  the  two-dimensipnal  theory  it  is  assumed  that 
the  energy  of  the  vortex,  E^,  is  equal  to  the  entrained  energy  E^  If  the  vortex  is  inclined  to 
the  oncoming  stream,  as  sketched  in  Figure  7,  then  not  all  of  the  entrained  energy  is 
converted  into  mixing  since  there  will  now  be  an  energy  flux  in  the  "spanwise”  direction 
through  the  mixing  layer.  Hence 


Ev  =  i|E,  (58) 

where  f)  represents  the  fraction  of  entrained  energy  traitsformed  into  rotational  energy. 
Combiiting  Eq.  (53),  (54),  and  (58)  and  solving  for  m  ^ves 


m 


y(y-i) 

2 


O  O 


(59) 


As  in  Ref.  2  Eq  Is  nondimensionalized  by  Ej^o^/  where 


E 

m 


(60) 


to  give 

i  =  Eo(Ej>oV) 

where  V  is  the  slowly  growing  control  volume  ABCD  in  Figure  6.  If  UgS  in  Eq.  (59)  is 
chosen  such  that 


U^t«V 


Then  Eq.  (59)  gives 


m 


//  ~ 
E. 


(61) 


(62) 


The  vortex  model  is  the  flow  over  a  two-dimensional  vortex  by  a  freestream  normal  to 
the  vortex.  In  Jhe  frame  of  reference  moving  with  the  vortex,  the  mixed  layer  has  a 
temperature  of  T  and  the  oncoming  normal  flow  has  a  Mach  number  of  M^os9  (see  Figure 
7).  Consequently,  the  reference  values,  T.  and  M.  are  chosen  to  be  as  foEows: 


T_  =  T 

M.=  M^os6 


(63) 
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q.  =  qjCOse  J 

where  <1.  is  the  convective  vdocity  of  the  mixing  layer  and  Eq.  (29)  now  becomes 

n  B  1}  — 


®o  [1  ♦ 


(64) 


It  now  remains  to  estimate  the  conversion  factor,  q.  If 

l^<<2/y-l 

theit,  to  a  first  approximation  (see  Eq.  (57)) 


(y~i) 

2y 


(65) 


and  consequently 


q  =  Ey/Eg  =  cos^6 


(66) 


Hence  the  final  result  is 


From  Eq.  (56)  and  (57) 


(67) 


K  *  1  ‘®®* 

V 

who-e  Aq  is  the  nondimensionalized  (by  qj  velocity  perturbation  due  to  the  vortex  and  V  is 
the  nondimensional  volume  V.  Hence, 


E 

O 


RKE 


(69) 


where  RKE  is  the  nondimensionalized  rotational  kinetic^eigy  of  the  vortex.  The  mixing 
can  therefore  be  found  in  terms  of  the  reference  mixing  m^  at^  the  reference  RKE,  (RKEIg. 
If  the  mixing  is  divided  by  its  value  at  denoted  by  the  subscript  i  (the  usual  method 

of  presenting  mixing  results),  then 
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The  first  term  in  braces  is  the  result  found  for  two-dimensional  vortices  in  Ref.  2,  and  by 
using  the  Prandtl-Glauert  equation,  Eq.  (SO),  it  was  shown  that 
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(RKB)oi 
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••/2 

(71) 


Using  the  definition  of  M.  from  Eq.  (63)  then 


(72) 


Eq.  (72)  represents  the  mixing  rate  variation  with  and  the  unlcnown  sweep  angle,  0.  It 

remains  now  to  determine  the  sweep  ang^e 

Simple  difierentiation  of  Eq.  (72)  shows  that  there  are  mbdi^  maximum  points  when 


=  0 

It  is  suggested  that  these  maximum  points  are  the  only  situations  that  can  exist  in  practice. 
They  occur,  for  y  =  1.4,  at 

e=o 

M^cosO  =  |0.662 
11.777 

For  Mf  <  .662  only  the  unswept  vortex  is  in  equilibrium.  If  a  .662,  the  vortex  angle,  0,  for 

maximum  mixing  is  given  by 

cosO  =  .662/Mc  (74) 

which  gives  two  values  of  0,  namely 

e  =  ±jcos'*(.662/M^|  (75) 
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The  third  possibility,  MeCos0  «  1.777,  is  not  of  present  interest  since  the  lowest  possible 
value  of  Mci^lsatisfies  this  relation  (M«  =  1.7^  does  not  remotely  satisfy  Eq.  (55).  The 
variation  of  m/nri)  with  9  for  various  Mach  numbers  is  shown  in  Rgure  8. 

The  question  now  arises  as  to  when  the  swept  modes  actually  appear  in  an  experiment 
It  is  platrsible  to  argue  that  a  smooth  transition  from  the  unswept  to  the  swept  mode  will 
only  occur  when  the  mixing  rate  is  the  same  for  both  modes,  since  this  would  not  involve  a 
discontinuous  energy  entrainment  Consequently,  it  is  suggested  that  until  a  value  of 
only  the  unswept  mode  is  apparent,  where  M^  is  the  value  for  which 


% 'unswept 


(76) 


and 


swept 


is  the  fflaxlmum  value 


given  by  (for  y  =  1.4) 


swept 


.21047 


1-  .05893 


(77) 


Eq.  (76)  is  satisfied  when 

M;=.662 

Eq.  (77)  is  shown  as  the  dashed  curve  in  Figure  2.  The  combined  two-dimensional  and 
swept  vortex  analysis  agree  qualitatively  with  observed  mixing  decreases  over  a  wide  range 
of  M^.  Its  quantitative  predictions  of  mixing  decrease  run  right  through  the  center  of  the 
data  scatter  almost  until  =  1,  after  which  the  theory  seems  to  slightly  underpredict  the 

observed  rates.  Note  that  the  swept  vortex  analysis  predicts  a  finite  mixing  rate  for  large 
and  that  this  rate  shoiild  be  reached  by  =  3.  Thus,  die  invocation  of  a  "mixing  maximum” 

principle  within  the  context  of  previously  described  energy  considerations  accounts 
qualitatively  and,  to  a  large  extent,  quantitatively  for  observed  mixing  decreases.  The 
theory  implies  that  the  flattening  of  mixing  rates  for  >  .66  should  be  associated  with  the 

appearance  of  swept  or  symmetrically  swept  ribbed  structures  in  the  flow. 


ROUND  JETS 

Once  again  the  mixing  is  assumed  dominated  by  vortices.  There  are  two  types  of 
observed  vortical  structures  in  the  round  jet,  namely  axisymmetric  and  helical  vortices.  In 
both  cases  the  mixing  is  given  by  the  integial 
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1 

m  =  p  I  v^dx  (78) 

-1 

where  v,  is  the  radial  velocity,  x  is  the  axial  direction  and  21  denotes  a  characteristic  length 
between  vortical  structures,  as  shown  in  Hgiue  9. 

Axisymmetric  Vortices 

The  velocity  induced  by  an  axisymmetric  vortex  ring,  with  radius  r^  is  governed  by  the 
Prandd-Glauert  equation 

r 

where  6  is  the  azimuthal  angle,  #  is  the  perturbation  velocity  potential,  is  the  convective 

Mach  number,  and 


p2  =  l-lV^ 


(80) 


For  axisymmetric  structures  die  azimuthal  variation  is  zero.  Using  Eq.  (79)  for  a  ring  vortex 
gives  the  following  formulae  for  the  axial  induced  velocity,  v^,  and  the  radial  induced 
velocity,  v^ 


V  s  _ - — 

r  4nr^ 


2X 


I 


2x 


1 


Cl  -  r  co8(e>P')]d0' 

TP - 3 -  - 

+  r  +  1  -  2r  co8(P-P')] 

X  co8(P-P'))de» 

- - 2 - - 

+  r  +  1  -  2r  co8{8-8')] 


(81) 


(82) 


where  r  is  the  strength  of  the  vortex  and  r  is  the  radius  nondimensionalized  by  r,,. 


Helical  Vortices 


The  effective  Mach  number  is  normal  to  the  vortex  and  hence  if  the  angle  of  die  hdix  is 
9,  the  eff^ve  Mach  number  is  M,.  cos  9,  and  therefore 

p2  =  i-M2cos29,  (83) 

The  induced  velocities  given  by  one  rotation  of  the  helix  (9  0, 2x)  are  as  follows: 
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wl^re  Vg  is  the  azimuthal  velocity  and 

ft  =  tantf,  (87) 

In  the  analysis  above  the  angle  of  the  helix  is  unknown,  but  fortunately  some 
experimental  data  is  available  from  the  acoustics  studies  reported  by  Seiner.^^  The  data 
makes  clear  that  several  wavelei^;ths  or  modes  can  exist  simultaneoudy  at  certain  values  of 
M|.  The  jet  Mach  number,  can  be  crudely  reiated  to  the  convective  Mach  number  M^  by 

M,*2Mc  (88) 


andthehdicalai^e,  9^  can  be  rdated  to  the  wave  length  by 


u  3  tan  9 

8 


A 

2nro 


(89) 


where  A  is  the  wavelength.  Seiner  reports  that  above  a  jet  Mach  number  of  1.2  the  vortical 
structures  are  helical.  It  is  assumed  that  the  dominant  source  of  noise  is  the  vortices  passing 
a  point  in  space 

The  variation  of  ft  with  M,.  can  be  found  from  a  simple  linear  fit  to  the  data  in  Ref.  13 
and  by  using  Equations  (88)  and  (89).  Thus 


H  = 


i  (1  +  44MJ  ^ow  mode" 

It 

~  (1.5  +  5.5MJ  "intermediate  mode" 

It 

I  (1.75  +  62MP  "high  mode" 


(90) 


One  of  the  most  interesting  features  of  these  angles  is  that 
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This  implies  tihat  the  Mach  number  normal  to  the  vortical  structure/  M^osO,  is  less  than 
0.71.  In  other  words,  the  oncoming  flow  is  subsoidc.  In  the  two-dimensional  studies 
reported  in  Reference  2  shock  waves  occurred  for  Mach  numbers  above  0.78,  implying  that 
for  the  rotmd  jet  the  structure  switches  from  axisymmetric  vortices  to  heUcal  vortices  to 
avoid  shock  waves,  a  hypothesis  first  advanced  in  Ref.  2. 

Mixing  ralmlarinm 

The  formulae  for  the  velocities  given  above  can  now  be  used  to  compute  the  mixing 
ratio  given  by  Eq.  (3).  In  these  calculations  the  effect  of  a  large  number  of  axisymmetric 
vortices  or  a  large  number  of  turns  of  the  helical  vortices  are  represented.  There  is  a 
problem  with  scaling  the  mixing  rate,  because  the  value  of  the  incompressible  mixing  is  not 
clear  from  the  Schadow  et  al*^  data.  Cot\sequently,  for  comparison  purposes,  the  following 
assumptions  are  made. 

(a)  For  Me  <  0.6  axisymmetric  vortices  are  assumed  to  exist  and  the  scaling  of 
the  mixing  is  fixed  to  agree  with  Schadow  et  al^*  data  at  M^  =  0.26. 

(b)  For  M^  >  0.6  helical  vortices  are  assumed  to  exist.  This  figure  is  obtained 
from  information  given  Seiner.^  The  scaling  is  fixed  to  match  the  data  at 

Me  =  0.6. 


Results  from  the  computations  are  shown  in  Figure  10  and  it  may  be  seen  that  theory 
agrees  surprisingly  well  with  the  data.  Also  ^own  are  two  curves  computed  using  the  high 
and  low  modes  for  the  helical  vortices:  these  are  not  tied  to  any  experimental  data  points 
and  are  scaled  by  the  value  a  helical  vortex  model  would  give  at  M,.  =  0.0.  It  can  be  seen 
that  mode  switching  produces  little  change  in  tlw  mixing.  The  most  promising  avenue  to 
increase  mixing  is  to  switch  from  an  axisymmetric  structure  to  a  hdical  structure  at  as  low  a 
Mach  number  as  possible. 


VARIATIONAL  APPROACH  TO  HEAT  RELEASE  MAXIMIZATION 
The  desire  is  to  maximize  the  volumetric  l^t  release  H 
Xo+L  - 

H  -  j  I  hj(#,  (92) 

x=Xo  y=— 

subject  to  the  constraint  that  the  velocity  field  satisfy  Eq.  (41).  The  key  here  is  to  find 
plausible  models  to  link  the  specific  heat  rdease  hf  to  the  local  velocity  field  (4, #y— )• 
Functionals  have  been  defined  dependent  upon  streamline  length  and/or  some  positive 
definite  function  of  the  streamline  curvature.  These  two  quahtities  provide  a 
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characterization  of  the  length  and  wrinkledness  of  interfaces  between  reactants,  and  the 
expectation  was  that  the  longer  and  more  convoluted  this  interfeice,  the  greater  the  total 
product  formation  and  heat  release.  The  variational  problem  involves  finding  a  flow 
configuration  that  balances  streamline  convolution  against  the  diffusive  or  averaging 
character  of  the  potential  flow  constraint 


The  first  problem  to  be  dealt  with  is  the  description  of  the  streamline  length  and 
curvature  in  terms  of  the  potential  #.  By  writing  the  formal  expressions  for  arclength  and 
curvature  of  the  streamfunction  |r(x,y)  in  terms  of  u  and  v  and  then  substituting  the 
appropriate  expressions  in  terms  of  the  potoitial  #,  one  obtains  the  length  element  as 

S  »  [  1  +  (  #y  /  #3,  (93) 

and  the  curvature  K  as 


K  = 


^  ♦x  [  ^  ^  ] 

[  1  ♦  (  ♦y  /  ♦x**  ] 


(94) 


Since  the  streamline  length  is  a  positive  definite  quantity  a  possible  relation  between  it  and 
heat  release  can  be  plausibly  assumed.  The  curvature,  on  the  other  hand,  may  be  negative, 
and  thus  the  heat  release  can  be  coimected  to  it  only  flmmgh  some  positive  valued  fimction 
of  K.  The  nmgnitude  of  K  is  an  obvious  choice,  but  leads  to  analytical  difficulties  in 
obtaining  the  appropriate  Euler  equation  for  the  variational  problem.  Instead,  it  has  been 
assumed  hf  -  »  F.  The  Euler  equation  for  this  functional  is 


d_ 

dx 


=  0 


(95) 


Expanding  the  x-  and  y-  partial  derivatives  gives 
’^xx  ’^xy  ’^yy  ’^xx’^xy 

u  xx-'w  »xy  j  »y»yy#yy^  ♦xyV 


r  ^  4 


•t-  ^  ^ 

xy  xxy  XX  xyy 


X  xy  xy  xy  xy  xy 


[44  +44  If 

»-  xy  xyy  yy  xxy  ■»  > 


^y^xy^xy 


=  0. 


(96) 


Realizing  that  all  of  the  partial  derivatives  of  F  with  respect  to  ^xy  rational 

functions  of  these  arguments,  it  is  quickly  clear  that  the  resulting  equation  for  4  is  high 
order  and  strongly  nonlinear(  6'*'  order  polynomials  in  4xf  obvious  simplifications 

have  been  found  to  make  the  equation  analytically  tractable,  and  the  strong  nonlinearity 
urges  caution  in  ntuneiical  treatment 
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Figure  1.  Vorticity  contours  in  a  two-dimensional  mixing  layer,  M|  =  2.0,  M2  =  1.2,  M, 


0*0 


Figure  3.  Spreading  rates  for  unheated  and  heated  two-dimensional  shear  layers,  M, 


Figure  4.  Normalized  temperature  contours  in  the  heated  two-dimensional  shear  layer,  M, 


Figure  5.  Entrainment  velocities  due  to  heat  release,  f  *  <  0. 


Figure  6.  Control  volume  for  energy  balance  in  a  shear  layer. 


Figure  ?.  Swept  vortices 


versus  9  for  various  convective  Mach  numbers,  M, 


